Single-Cell Transcriptional Profile Construction of Rat Pituitary Glands before and after Sexual Maturation and Identification of Novel Marker Spp1 in Gonadotropes.

The mammalian pituitary gland drives highly conserved physiological processes such as somatic cell growth, pubertal transformation, fertility, and metabolism by secreting a variety of hormones. Recently, single-cell transcriptomics techniques have been used in pituitary gland research. However, more studies have focused on adult pituitary gland tissues from different species or different sexes, and no research has yet resolved cellular differences in pituitary gland tissue before and after sexual maturation. Here, we identified a total of 15 cell clusters and constructed single-cell transcriptional profiles of rats before and after sexual maturation. Furthermore, focusing on the gonadotrope cluster, 106 genes were found to be differentially expressed before and after sexual maturation. It was verified that Spp1, which is specifically expressed in gonadotrope cells, could serve as a novel marker for this cell cluster and has a promotional effect on the synthesis and secretion of follicle-stimulating hormone. The results provide a new resource for further resolving the regulatory mechanism of pituitary gland development and pituitary hormone synthesis and secretion.

Characterization of sexual maturity-associated N6-methyladenosine in boar testes.

BACKGROUND: The health and size of the testes are crucial for boar fertility. Testicular development is tightly regulated by epigenetics. N6-methyladenosine (m6A) modification is a prevalent internal modification on mRNA and plays an important role in development. The mRNA m6A methylation in boar testicular development still needs to be investigated. RESULTS: Using the MeRIP-seq technique, we identify and profile m6A modification in boar testes between piglets and adults. The results showed 7783 distinct m6A peaks in piglets and 6590 distinct m6A peaks in adults, with 2,471 peaks shared between the two groups. Enrichment of GO and KEGG analysis reveal dynamic m6A methylation in various biological processes and signalling pathways. Meanwhile, we conjointly analyzed differentially methylated and expressed genes in boar testes before and after sexual maturity, and reproductive related genes (TLE4, TSSK3, TSSK6, C11ORF94, PATZ1, PHLPP1 and PAQR7) were identified. Functional enrichment analysis showed that differential genes are associated with important biological functions, including regulation of growth and development, regulation of metabolic processes and protein catabolic processes. CONCLUSION: The results demonstrate that m6A methylation, differential expression and the related signalling pathways are crucial for boar testicular development. These results suggest a role for m6A modification in boar testicular development and provided a resource for future studies on m6A function in boar testicular development.

Loss of Fshr Prevents Testicular Maturation in Atlantic Salmon (Salmo salar L.).

Early puberty poses a significant challenge for male Atlantic salmon in aquaculture due to its negative impact on growth and welfare. The regulation of puberty in vertebrates involves 2 key reproductive hormones: follicle-stimulating hormone (FSH) and luteinizing hormone (LH) and their gonadal receptors. In male mice lacking FSH receptor, testes size is reduced, but fertility is maintained, while medaka and zebrafish with a disrupted fshr gene exhibit near normal testis size and fertility. In these fishes both Fsh and Lh are present during puberty and Lh may rescue fertility, while in salmonid fish only Fsh is present in the circulation during puberty. Using CRISPR-Cas9, we produced crispants with a high prevalence of fshr mutations at the target site, which remained fertile, although more than half showed a testis development deviating from wild-type (wt) males. Crossing out these F0 crispants to each other produced a viable F1 generation showing frameshift (fshr-/-) or in-frame mutations (fshrif/if). Nearly all wt males matured while all fshr-/- males remained immature with small testes containing A spermatogonia as the furthest developed germ cell type and prepubertal plasma androgen levels. Also, the pituitary transcript levels of gnrhr2bba and lhb, but not for fshb, were reduced in the fshr-/- males compared with maturing males. More than half of the fshrif/if mutant males showed no or a delayed maturation. In conclusion, Atlantic salmon show the unique characteristic that loss of Fshr function alone results in male infertility, offering new opportunities to control precocious puberty or fertility in salmon.

Probabilistic reaction norm reveals family-related variation in the association between size, condition, and sexual maturation onset in Atlantic salmon (Salmo salar).

This study investigated the relationship between the size, condition, year class, family, and sexual maturity of Atlantic salmon (Salmo salar) using data collected in an aquaculture selective breeding programme. Males that were sexually mature at 2 years of age (maiden spawn) have, on average, greater fork length and condition factor (K) at 1 year of age than their immature counterparts. For every 10-mm increase in fork length or 0.1 increase in K at 1 year of age, the odds of sexual maturity at 2 years of age increased by 1.48 or 1.22 times, respectively. Females that were sexually mature at 3 years of age (maiden spawn) have, on average, greater fork length and K at 2 years of age than their immature counterparts. For every 10-mm increase in fork length or 0.1 increase in K at 2 years of age, the odds of sexual maturity at 3 years of age increased by 1.06 or 1.44 times, respectively. The family explained 34.93% of the variation in sexual maturity among 2-year-old males that was not attributable to the average effects of fork length and K at 1 year of age and year class. The proportion of variation in sexual maturity among 3-year-old females explained by the family could not be investigated. These findings suggest that the onset of sexual maturation in Atlantic salmon is conditional on performance (with respect to energy availability) surpassing a threshold, the magnitude of which can vary between families and is determined by a genetic component. This could support the application of genetic selection to promote or inhibit the onset of sexual maturation in farmed stocks.

Clustering of multi-tissue transcriptomes in gilts with normal cyclicity or delayed puberty reveals genes related to pubertal development†.

In gilts, puberty is marked by standing estrus in the presence of a boar. Delayed puberty (DP; failure to display pubertal estrus) is a major reason for gilt removal. To investigate the physiological determinants underlying DP in gilts, transcriptomic data from tissues relevant to estrus and puberty, such as mediobasal hypothalamus, anterior pituitary gland, ovarian cortex, olfactory bulb, amygdala, and hippocampus, were obtained from age-matched DP (n = 8) and cyclic control gilts at follicular phase (n = 8) and luteal phase (n = 8) of the estrous cycle. A gene expression module analysis via three-way gene × individual × tissue clustering using tensor decomposition identified pituitary and ovary gene modules contributing to regulation of pubertal development. Analysis of gene expression in the hypothalamic-pituitary-ovary axis identified reduced expression of hypothalamic genes critical for stimulating gonadotropin secretion (KISS1 and TAC3) and reduced expression of LHB in the anterior pituitary of DP gilts compared with their cyclic counterparts. Consequently, luteinizing hormone-induced genes in the ovary important for folliculogenesis (OXTR, RUNX2, and PTX3) were less expressed in DP gilts. Other intrafollicular genes (AHR, PTGS2, PTGFR, and IGFBP7) and genes in the steroidogenesis pathways (STAR and CYP11A1) necessary to complete the ovulatory cascade were also less expressed in DP gilts. This is the first clustering of multi-tissue expression data from DP and cyclic gilts to identify genes differentially expressed in gilts of similar ages but at different levels of sexual development. A critical lack of gonadotropin support and reduced ovarian responsiveness underlie DP in gilts.

Regulatory Effects of the Kiss1 Gene in the Testis on Puberty and Reproduction in Hezuo and Landrance Boars.

Kisspeptin, a neuropeptide encoded by the Kiss1 gene, combines with its receptor Kiss1R to regulate the onset of puberty and male fertility by the hypothalamic-pituitary-gonadal axis. However, little is known regarding the expression signatures and molecular functions of Kiss1 in the testis. H&E staining revealed that well-arranged spermatogonia, spermatocytes, round and elongated spermatids, and spermatozoa, were observed in 4-, 6-, and 8-month-old testes compared to 1- and 3-month-old testes of Hezuo pigs; however, these were not observed in Landrance until 6 months. The diameter, perimeter, and cross-sectional area of seminiferous tubules and the perimeter and area of the tubular lumen increased gradually with age in both pigs. Still, Hezuo pigs grew faster than Landrance. The cloning results suggested that the Hezuo pigs' Kiss1 CDS region is 417 bp in length, encodes 138 amino acids, and is highly conserved in the kisspeptin-10 region. qRT-PCR and Western blot indicated that the expression trends of Kiss1 mRNA and protein were essentially identical, with higher expression levels at post-pubertal stages. Immunohistochemistry demonstrated that the Kiss1 protein was mainly located in Leydig cells and post-pubertal spermatogenic cells, ranging from round spermatids to spermatozoa. These studies suggest that Kiss1 is an essential regulator in the onset of puberty and spermatogenesis of boars.

Molecular characterization of Smtdc-1 and Smddc-1 discloses roles as male-competence factors for the sexual maturation of Schistosoma mansoni females.

Introduction: Schistosomes are the only mammalian flatworms that have evolved separate sexes. A key question of schistosome research is the male-dependent sexual maturation of the female since a constant pairing contact with a male is required for the onset of gonad development in the female. Although this phenomenon is long known, only recently a first peptide-based pheromone of males was identified that contributes to the control of female sexual development. Beyond this, our understanding of the molecular principles inducing the substantial developmental changes in a paired female is still rudimentary. Objectives: Previous transcriptomic studies have consistently pointed to neuronal genes being differentially expressed and upregulated in paired males. These genes included Smp_135230 and Smp_171580, both annotated as aromatic-L-amino-acid decarboxylases (DOPA decarboxylases). Here, we characterized both genes and investigated their roles in male-female interaction of S. mansoni. Methodologies/findings: Sequence analyses indicated that Smp_135230 represents an L-tyrosine decarboxylase (Smtdc-1), whereas Smp_171580 represents a DOPA decarboxylase (Smddc-1). By qRT-PCR, we confirmed the male-specific and pairing-dependent expression of both genes with a significant bias toward paired males. RNA-interference experiments showed a strong influence of each gene on gonad differentiation in paired females, which was enhanced by double knockdown. Accordingly, egg production was significantly reduced. By confocal laser scanning microscopy, a failure of oocyte maturation was found in paired knockdown females. Whole-mount in situ hybridization patterns exhibited the tissue-specific occurrence of both genes in particular cells at the ventral surface of the male, the gynecophoral canal, which represents the physical interface of both genders. These cells probably belong to the predicted neuronal cluster 2 of S. mansoni. Conclusion: Our results suggest that Smtdc-1 and Smddc-2 are male-competence factors that are expressed in neuronal cells at the contact zone between the genders as a response of pairing to subsequently control processes of female sexual maturation.

Genotype-specific variation in seasonal body condition at a large-effect maturation locus.

Organisms use resource allocation strategies to survive seasonal environmental changes and life-history stage transitions. Earlier studies found a transcription cofactor, vgll3, associating with maturation timing that inhibits adipogenesis in mice and affects body condition in juvenile salmon. Owing to a lack of temporal studies examining seasonality effects on phenotypes such as vgll3 genotype, body condition, maturation and different life stages, we investigated the influence of different larval and juvenile temperatures, vgll3 genotype and interactions with body condition and maturation rate. We reared Atlantic salmon for 2 years in four larval-juvenile phase temperature groups until the occurrence of mature males. We found no effect of larval temperature on the measured phenotypes or maturation rate. However, we observed an increased maturation rate in individuals of the warm juvenile temperature treatment and differences in body condition associated with vgll3 genotype. Early maturation genotype individuals had a less variable body condition across seasons compared with late maturation genotype individuals. This result suggests a vgll3 influence on resource allocation strategies; possibly linked with the early maturation process, with early maturation genotype individuals having a higher maturation rate and a higher body condition in the spring.

Enduring sex-dependent effects of lipopolysaccharide treatment on the hypothalamic-pituitary-gonadal axis in mice.

Pubertal stress causes enduring sexual behavior dysfunction in males and females, but the underlying mechanism remains unknown. These changes may arise from pubertal programming of the hypothalamic-pituitary-gonadal axis. Previous findings show that stress exposure downregulates the hypothalamic-pituitary-gonadal axis, particularly through the reduction of the neuropeptide kisspeptin (Kiss1) and its receptor (Kiss1R). Although acute changes in kiss1 and Kiss1r genes have been observed following pubertal immune stress, it is unclear whether immune stress-induced downregulation of kiss1 and kiss1r persists beyond puberty. The current study investigated the enduring sex-specific consequences of lipopolysaccharide on the expression of Kiss1 and Kiss1r in 160 pubertal or adult mice at multiple time points. Six-week and 10-week-old male and female mice were treated with either saline or with lipopolysaccharide. Mice were euthanized either 8 h or 4 weeks following treatment. Although we did not identify any sex differences, our results revealed that lipopolysaccharide treatment decreases hypothalamic Kiss1 and Kiss1r in both pubertal and adult mice within 8 h of treatment. The decreased hypothalamic Kiss1 expression persists 4 weeks later only in mice treated with lipopolysaccharide during puberty. Our findings highlight the age-dependent vulnerability of the hypothalamic-pituitary-gonadal axis to immune stress, providing a better understanding of the mechanisms implicated in allostatic shift during immune stress. Finally, our findings also show the effects of immune stress on various components of the hypothalamic-pituitary-gonadal axis, which could have implications for sexual and fertility-related dysfunctions.

gBLUP-GWAS identifies candidate genes, signaling pathways, and putative functional polymorphisms for age at puberty in gilts.

Successful development of replacement gilts determines their reproductive longevity and lifetime productivity. Selection for reproductive longevity is challenging due to low heritability and expression late in life. In pigs, age at puberty is the earliest known indicator for reproductive longevity and gilts that reach puberty earlier have a greater probability of producing more lifetime litters. Failure of gilts to reach puberty and display a pubertal estrus is a major reason for early removal of replacement gilts. To identify genomic sources of variation in age at puberty for improving genetic selection for early age at puberty and related traits, gilts (n = 4,986) from a multigeneration population representing commercially available maternal genetic lines were used for a genomic best linear unbiased prediction-based genome-wide association. Twenty-one genome-wide significant single nucleotide polymorphisms (SNP) located on Sus scrofa chromosomes (SSC) 1, 2, 9, and 14 were identified with additive effects ranging from -1.61 to 1.92 d (P < 0.0001 to 0.0671). Novel candidate genes and signaling pathways were identified for age at puberty. The locus on SSC9 (83.7 to 86.7 Mb) was characterized by long range linkage disequilibrium and harbors the AHR transcription factor gene. A second candidate gene on SSC2 (82.7 Mb), ANKRA2, is a corepressor for AHR, suggesting a possible involvement of AHR signaling in regulating pubertal onset in pigs. Putative functional SNP associated with age at puberty in the AHR and ANKRA2 genes were identified. Combined analysis of these SNP showed that an increase in the number of favorable alleles reduced pubertal age by 5.84 ± 1.65 d (P < 0.001). Candidate genes for age at puberty showed pleiotropic effects with other fertility functions such as gonadotropin secretion (FOXD1), follicular development (BMP4), pregnancy (LIF), and litter size (MEF2C). Several candidate genes and signaling pathways identified in this study play a physiological role in the hypothalamic-pituitary-gonadal axis and mechanisms permitting puberty onset. Variants located in or near these genes require further characterization to identify their impact on pubertal onset in gilts. Because age at puberty is an indicator of future reproductive success, these SNP are expected to improve genomic predictions for component traits of sow fertility and lifetime productivity expressed later in life.

Selecting for replacement gilts is challenging because sow reproductive traits are lowly heritable and expressed late in life. Age at puberty is the earliest indicator of future reproductive success of gilts. Genetic selection for early onset of puberty could be feasible with the availability of molecular genetic predictors for age at puberty. To identify genomic sources associated with variation in age at puberty in gilts, a large-scale genome-wide association study was conducted at the U.S Meat Animal Research Center, Clay Center, Nebraska. Novel genomic associations for age at puberty were identified. Several candidate genes identified for age at puberty are involved in signaling pathways that regulate ovarian functions and pubertal onset. Potential causative genetic variants for age at puberty were identified within the candidate genes. These novel SNP are important new markers for use in genomic selection of replacement gilts with early puberty and provide critical new insight into biological mechanisms important for pubertal development in gilts.

A pituitary gene network linking vgll3 to regulators of sexual maturation in male Atlantic salmon.

Age at maturity is a key life history trait and a significant contributor to life history strategy variation. The maturation process is influenced by genetic and environmental factors, but specific causes of variation in maturation timing remain elusive. In many species, the increase in the regulatory gonadotropin-releasing hormone 1 (GnRH1) marks the onset of puberty. Atlantic salmon, however, lacks the gnrh1 gene, suggesting gnrh3 and/or other regulatory factors are involved in the maturation process. Earlier research in Atlantic salmon has found a strong association between alternative alleles of vgll3 and maturation timing. Recently we reported strong induction of gonadotropin genes (fshb and lhb) in the pituitary of Atlantic salmon homozygous for the early maturation allele (E) of vgll3. The induction of gonadotropins was accompanied by increased expression of their direct upstream regulators, c-jun and sf1 (nr5a1b) but the regulatory connection between vgll3 and these regulators has never been investigated in any organism. In this study, we investigated the potential regulatory connection between vgll3 genotypes and these regulators through a stepwise approach of identifying a gene regulatory network (GRN) containing c-jun and sf1, and transcription factor motif enrichment analysis. We found a GRN containing c-jun with predicted upstream regulators, e2f1, egr1, foxj1 and klf4, to be differentially expressed in the pituitary. Finally, we suggest a vgll3 and Hippo pathway -dependent model for transcriptional regulation of c-jun and sf1 in the pituitary, which may have broader implications across vertebrates.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits.

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

Genetic Basis of Sexual Maturation Heterosis: Insights From Ovary lncRNA and mRNA Repertoire in Chicken.

Sexual maturation is fundamental to the reproduction and production performance, heterosis of which has been widely used in animal crossbreeding. However, the underlying mechanism have long remained elusive, despite its profound biological and agricultural significance. In the current study, the reciprocal crossing between White Leghorns and Beijing You chickens were performed to measure the sexual maturation heterosis, and the ovary lncRNAs and mRNAs of purebreds and crossbreeds were profiled to illustrate molecular mechanism of heterosis. Heterosis larger than 20% was found for pubic space and oviduct length, whereas age at first egg showed negative heterosis in both crossbreeds. We identified 1170 known lncRNAs and 1994 putative lncRNAs in chicken ovary using a stringent pipeline. Gene expression pattern showed that nonadditivity was predominant, and the proportion of nonadditive lncRNAs and genes was similar between two crossbreeds, ranging from 44.24% to 49.15%. A total of 200 lncRNAs and 682 genes were shared by two crossbreeds, respectively. GO and KEGG analysis showed that the common genes were significantly enriched in the cell cycle, animal organ development, gonad development, ECM-receptor interaction, calcium signaling pathway and GnRH signaling pathway. Weighted gene co-expression network analysis (WGCNA) identified that 7 out of 20 co-expressed lncRNA-mRNA modules significantly correlated with oviduct length and pubic space. Interestingly, genes harbored in seven modules were also enriched in the similar biological process and pathways, in which nonadditive lncRNAs, such as MSTRG.17017.1 and MSTRG.6475.20, were strongly associated with nonadditive genes, such as CACNA1C and TGFB1 to affect gonad development and GnRH signaling pathway, respectively. Moreover, the results of real-time quantitative PCR (RT-qPCR) correlated well with the transcriptome data. Integrated with positive heterosis of serum GnRH and melatonin content detected in crossbreeds, we speculated that nonadditive genes involved in the GnRH signaling pathway elevated the gonad development, leading to the sexual maturation heterosis. We characterized a systematic landscape of ovary lncRNAs and mRNAs related to sexual maturation heterosis in chicken. The quantitative exploration of hybrid transcriptome changes lays foundation for genetic improvement of sexual maturation traits and provides insights into endocrine control of sexual maturation.

Dicer ablation in Kiss1 neurons impairs puberty and fertility preferentially in female mice.

Kiss1 neurons, producing kisspeptins, are essential for puberty and fertility, but their molecular regulatory mechanisms remain unfolded. Here, we report that congenital ablation of the microRNA-synthesizing enzyme, Dicer, in Kiss1 cells, causes late-onset hypogonadotropic hypogonadism in both sexes, but is compatible with pubertal initiation and preserved Kiss1 neuronal populations at the infantile/juvenile period. Yet, failure to complete puberty and attain fertility is observed only in females. Kiss1-specific ablation of Dicer evokes disparate changes of Kiss1-cell numbers and Kiss1/kisspeptin expression between hypothalamic subpopulations during the pubertal-transition, with a predominant decline in arcuate-nucleus Kiss1 levels, linked to enhanced expression of its repressors, Mkrn3, Cbx7 and Eap1. Our data unveil that miRNA-biosynthesis in Kiss1 neurons is essential for pubertal completion and fertility, especially in females, but dispensable for initial reproductive maturation and neuronal survival in both sexes. Our results disclose a predominant miRNA-mediated inhibitory program of repressive signals that is key for precise regulation of Kiss1 expression and, thereby, reproductive function.

Hypothalamic Overexpression of Makorin Ring Finger Protein 3 Results in Delayed Puberty in Female Mice.

Makorin ring finger protein 3 (MKRN3) is an important neuroendocrine player in the control of pubertal timing and upstream inhibitor of gonadotropin-releasing hormone secretion. In mice, expression of Mkrn3 in the hypothalamic arcuate and anteroventral periventricular nucleus is high early in life and declines before the onset of puberty. Therefore, we aimed to explore if the persistence of hypothalamic Mkrn3 expression peripubertally would result in delayed puberty. Female mice that received neonatal bilateral intracerebroventricular injections of a recombinant adeno-associated virus expressing Mkrn3 had delayed vaginal opening and first estrus compared with animals injected with control virus. Subsequent estrous cycles and fertility were normal. Interestingly, male mice treated similarly did not exhibit delayed puberty onset. Kiss1, Tac2, and Pdyn mRNA levels were increased in the mediobasal hypothalamus in females at postnatal day 28, whereas kisspeptin and neurokinin B protein levels in the arcuate nucleus were decreased, following Mkrn3 overexpression, compared to controls. Cumulatively, these data suggest that Mkrn3 may directly or indirectly target neuropeptides of Kiss1 neurons to degradation pathways. This mouse model suggests that MKRN3 may be a potential contributor to delayed onset of puberty, in addition to its well-established roles in central precocious puberty and the timing of menarche.

Leveraging transcriptome and epigenome landscapes to infer regulatory networks during the onset of sexual maturation.

BACKGROUND: Despite sexual development being ubiquitous to vertebrates, the molecular mechanisms underpinning this fundamental transition remain largely undocumented in many organisms. We designed a time course experiment that successfully sampled the period when Atlantic salmon commence their trajectory towards sexual maturation. RESULTS: Through deep RNA sequencing, we discovered key genes and pathways associated with maturation in the pituitary-ovarian axis. Analyzing DNA methylomes revealed a bias towards hypermethylation in ovary that implicated maturation-related genes. Co-analysis of DNA methylome and gene expression changes revealed chromatin remodeling genes and key transcription factors were both significantly hypermethylated and upregulated in the ovary during the onset of maturation. We also observed changes in chromatin state landscapes that were strongly correlated with fundamental remodeling of gene expression in liver. Finally, a multiomic integrated analysis revealed regulatory networks and identified hub genes including TRIM25 gene (encoding the estrogen-responsive finger protein) as a putative key regulator in the pituitary that underwent a 60-fold change in connectivity during the transition to maturation. CONCLUSION: The study successfully documented transcriptome and epigenome changes that involved key genes and pathways acting in the pituitary - ovarian axis. Using a Systems Biology approach, we identified hub genes and their associated networks deemed crucial for onset of maturation. The results provide a comprehensive view of the spatiotemporal changes involved in a complex trait and opens the door to future efforts aiming to manipulate puberty in an economically important aquaculture species.

Strong regulatory effects of vgll3 genotype on reproductive axis gene expression in juvenile male Atlantic salmon.

Age at maturity is a major contributor to the diversity of life history strategies in organisms. The process of maturation is influenced by both genetics and the environment, and includes changes in levels of sex hormones and behavior, but the specific factors leading to variation in maturation timing are poorly understood. gnrh1 regulates the transcription of gonadotropin genes at pubertal onset in many species, but this gene is lacking in certain teleost species including Atlantic salmon (Salmo salar), which raises the possibility of the involvement of other important regulatory factors during this process. Earlier research has reported a strong association of alternative alleles of the vgll3 gene with maturation timing in Atlantic salmon, suggesting it as a potential candidate regulating reproductive axis genes. Here, we investigated the expression of reproductive axis genes in one-year-old Atlantic salmon males with immature gonads and different vgll3 genotypes during the spawning period. We detected strong vgll3 genotype-dependent differential expression of reproductive axis genes (such as fshb, lhb, amh and igf3) tested in the pituitary, and testis. In addition, we observed differential expression of jun (ap1) and nr5a1b (sf1), potential upstream regulators of gonadotropins in the pituitary, as well as axin2, id3, insl3, itch, ptgs2a and ptger4b, the downstream targets of amh and igf3 in the testis. Hereby, we provide evidence of strong vgll3 genotype-dependent transcriptional regulation of reproductive axis genes prior to sexual maturation and suggest alternative models for distinct actions of vgll3 genotypes on the related molecular processes.

Developmental sex differences in the peri-pubertal pattern of hypothalamic reproductive gene expression, including Kiss1 and Tac2, may contribute to sex differences in puberty onset.

The mechanisms regulating puberty still remain elusive, as do the underlying causes for sex differences in puberty onset (girls before boys) and pubertal disorders. Neuroendocrine puberty onset is signified by increased pulsatile GnRH secretion, yet how and when various upstream reproductive neural circuits change developmentally to govern this process is poorly understood. We previously reported day-by-day peri-pubertal increases (Kiss1, Tac2) or decreases (Rfrp) in hypothalamic gene expression of female mice, with several brain mRNA changes preceding external pubertal markers. However, similar pubertal measures in males were not previously reported. Here, to identify possible neural sex differences underlying sex differences in puberty onset, we analyzed peri-pubertal males and directly compared them with female littermates. Kiss1 expression in male mice increased over the peri-pubertal period in both the AVPV and ARC nuclei but with lower levels than in females at several ages. Likewise, Tac2 expression in the male ARC increased between juvenile and older peri-pubertal stages but with levels lower than females at most ages. By contrast, both DMN Rfrp expressionand Rfrp neuronal activation strongly decreased in males between juvenile and peri-pubertal stages, but with similar levels as females. Neither ARC KNDy neuronal activation nor Kiss1r expression in GnRH neurons differed between males and females or changed with age. These findings delineate several peri-pubertal changes in neural populations in developing males, with notable sex differences in kisspeptin and NKB neuron developmental patterns. Whether these peri-pubertal hypothalamic sex differences underlie sex differences in puberty onset deserves future investigation.

Pituitary Gonadotropin Gene Expression During Induced Onset of Postsmolt Maturation in Male Atlantic Salmon: In Vivo and Tissue Culture Studies.

Precocious male maturation causes reduced welfare and increased production costs in Atlantic salmon (Salmo salar) aquaculture. The pituitary produces and releases follicle-stimulating hormone (Fsh), the gonadotropin triggering puberty in male salmonids. However, little is known about how Fsh production is regulated in Atlantic salmon. We examined, in vivo and ex vivo, transcriptional changes of gonadotropin-related genes accompanying the initial steps of testis maturation, in pituitaries of males exposed to photoperiod and temperature conditions promoting maturation (constant light and 16°C). Pituitary fshb, lhb and gnrhr2bba transcripts increased in vivo in maturing males (gonado-somatic index > 0.1%). RNA sequencing (RNAseq) analysis using pituitaries from genetically similar males carrying the same genetic predisposition to mature, but differing by responding or not responding to stimulatory environmental conditions, revealed 144 differentially expressed genes, ~2/3rds being up-regulated in responders, including fshb and other pituitary hormones, steroid-related and other puberty-associated transcripts. Functional enrichment analyses confirmed gene involvement in hormone/steroid production and gonad development. In ex vivo studies, whole pituitaries were exposed to a selection of hormones and growth factors. Gonadotropin-releasing hormone (Gnrh), 17ß-estradiol (E2) and 11-ketotestosterone (11-KT) up-regulated gnrhr2bba and lhb, while fshb was up-regulated by Gnrh but down-regulated by 11-KT in pituitaries from immature males. Also pituitaries from maturing males responded to Gnrh and sex steroids by increased gnrhr2bba and lhb transcript levels, but fshb expression remained unchanged. Growth factors (inhibin A, activin A and insulin-like growth factor 1) did not change gnrhr2bba, lhb or fshb transcript levels in pituitaries either from immature or maturing males. Additional pituitary ex vivo studies on candidates identified by RNAseq showed that these transcripts were preferentially regulated by Gnrh and sex steroids, but not by growth factors, and that Gnrh/sex steroids were less effective when incubating pituitaries from maturing males. Our results suggest that a yet to be characterized mechanism up-regulating fshb expression in the salmon pituitary is activated in response to stimulatory environmental conditions prior to morphological signs of testis maturation, and that the transcriptional program associated with this mechanism becomes unresponsive or less responsive to most stimulators ex vivo once males had entered pubertal developmental in vivo.